Dark matter annihilation near a black hole: plateau vs. weak cusp 
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Dark matter annihilation in so-called "spikes" near black holes is believed to be an important method of 
indirect dark matter detection. In the case of circular particle orbits, the density profile of dark matter has 
a plateau at small radii, the maximal density being limited by the annihilation cross-section. However, in 
the general case of arbitrary velocity anisotropy the situation is different. Particulary, for isotropic velocity 
distribution the density profile cannot be shallower than r^^^^ in the very center. Indeed, a detailed study 
reveals that in many cases the term "annihilation plateau" is misleading, as the density actually continues to 
rise towards small radii and forms a weak cusp, p oc r~^^'^^^^\ where (3 is the anisotropy coefficient. The 
annihilation flux, however, does not change much in the latter case, if averaged over an area larger than the 
annihilation radius. 
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I. INTRODUCTION 



mass). Then Eq.([TJ becomes 



One of the most promising methods of indirect dark matter 
search is the detection of annihilation signal from so-called 
dark matter "spikes" around a black hole (either a supermasive 
black hole in galactic center or an intermediate-mass black 
hole). The existence of spikes is predicted by adiabatic growth 
model If the initial density profile is cusped, p (x , 
then the adiabatically grown spike will also have power-law 
profile p' oc , with 7' = (9 — 27) /(4 — 7). If the initial 
profile is cored, then the spike will have p' oc r^-^/"^. 

It is evident that densities can reach very high values for 
small radii (but greater than Vh ~ 2GMbh/ c^)- The annihi- 
lation rate of dark matter particles is therefore high enough 
to reduce the density. The usual argument is the following: 
consider the dark matter distribution function /(r,v). The 
annihilation rate at a certain point is given by 
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Here is the particle mass, 



is the annihilation cross- 



section times relative velocity, which is independent of v. 

If particles have circular orbits and the density is spherically 
symmetric, then one may integrate / over dv and obtain 
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Here po — p{r, < = 0) is the initial density, and pa is called 
annihilation plateau density. One can see that indeed for 
small radii the density approaches pa, and the "height" of the 
plateau decreases in time. 

In the general case of non-circular orbits (but still having 
spherical symmetry) it is easier to switch from {r, v} variables 
to {E, L} variables (energy and angular momentum per unit 
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p is the orbit-averaged density. For convenience, we replace 
L for i? = L'^/Ll, where Lc = GMth/V-^E is the angu- 
lar momentum of a circular orbit (and hence ^ i? ^ 1). 
The radius of a circular orbit is Vc — GMth/ {—2E), and we 
substitute r = xrr- Then we have 



p{xrc) 



dx 
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Assume we have density profile p{r) oc r '', and a velocity 
anisotropy parameter /3 — l~af/2a^ iQl- The case /3 = — 00 
corresponds to circular orbits, /3 = is the isotropic case, 
< /3 < 1 is the case of radial velocity anisotropy. 

The caveat is that in the isotropic case the density profile 
cannot be shallower than r~"^ in the center (or, more gener- 
ally, 7 ^ P + 1/2 id, Hi)- This means that we cannot get an- 
nihilation plateau with constant density for initially isotropic 
velocity distribution. One might argue that the annihilation in- 
troduces tangential velocity anisotropy to the degree compati- 
ble with flat density distribution, but in fact p(i?) is increasing 
function of R for < 7 < 1, and the circular orbits are in- 
deed depopulated. Thus, a careful examination is needed in 
the general case, as noted in |5]. This is the aim of the present 
study. 



II. A QUALITATIVE ARGUMENT 

A simple argument can explain the impossibility to create 
density profile more shallow than r^(>^+^/^'). 

Consider the distribution function of the following broken 
power-law form: 



fiE,R) = foR- 



{E/EoY' 
{E/E^Y- 



\E\ > \Eo\ 
\Eo\ > |£;| >0. 
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We assume that pi < p2, i.e. the distribution function is con- 
vex in log£' — logr coordinates. It is a simple exercise to 
show that the velocity anisotropy parameter equals exactly /3 
in the above expression. 

We now demonstrate that if pi > (3—1, then the den- 
sity in the region r tq = GMfjiJ (~Eq) is determined by 
the distribution function at \E\ > \Eo\, and p cx r'^^i+^Z^' 
(steeper than 7-^(^+1/2)). In the opposite case, however, the 
density in the region r < rp is determined by distribution 
function outside that region (i.e. with < |i?o|), provided 
thatp2 - /3 > -1. 

The density is given by expression 



fGMbhy^ 
dRf{ 
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Here e — —Er/GMbh is dimensionless energy. 

If we are interested in r <C ro, then the integral can be 
split in two constituents: pi — J^^ ■ ■ ■ and p2 — J^" • • • , 
representing the contribution of inner and outer areas, corre- 
spondingly. £0 = —Ecr/GMbii <C 1 by condition. Then we 
get 



Pi = 



3/2 



(7) 



1, then 

(3/2+pi)^ jjj 



X / de(i-e)-/5+i/2£P.-/5^-P. ^ 
i = 1,2, 

where limits of integration are given above. 

Now there are two cases for pi: if pi ~ (3 > 
the integrand is finite as Eq ^ and pi cx r 
the opposite case the integral diverges as Eq^^'^^^ and pi cx 
7.-(i/2+/3)^ The second integral, p2, is always cx 7-^(1/2+^). 
Hence in the first case pi ^ p2 and p is determined by pi, 
while in the second case they have the same dependence on r 
andpcxr-(^+i/2). 

What has this to do with annihilations? In the case /? < 
— 1 /2 we may get density plateau in the center due to annihila- 
tions, which corresponds to pi — —3/2. However, if initially 
(3 ^ —1/2, a constant density core cannot develop. Instead a 
sort of broken power-law density profile will emerge: 



r-(P2+3/2) 



r <ro 
r > ro 



(8) 



The break radius tq is the same as for density plateau, i.e. the 
density profile outside this radius is pa (r/ro)^*-'"^^'^/^^. Par- 
ticles outside ro are still not much affected by annihilation, 
since for them p ^ p{r) < pa- Inside ro, or, say another way, 
for > |i?o|, the distribution function is depleded more 
rapidly than in the case of annihilation core, since the aver- 
age density that particle "feels" is higher. The boundary ro 
corresponds to the intermediate area where p w Pa{t)- 
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FIG. 1 : Density profiles after a certain time of evolution ti (see text) 
for different values of /?: dotted line - circular orbits (/3 = —00); 
solid line ~ P — —2, dotted-dashed line - /3 = 0. In the latter case a 
weak cusp p cx r^^^^ develops instead of constant-density plateau. 
Initial density profile p oc r"^'^* is shown by the long-dashed line. 
Radius and density are scaled to the black hole influence radius and 
corresponding density. Annihilation plateau density is pa ~ 10®, 
corresponding annihilation radius ra = 2.2 ■ 10^'' is denoted by the 
vertical line. 



III. NUMERICAL CALCULATIONS 

The above qualitative arguments have to be confirmed by 
strict calculations. In order to do this, we solve the system of 
equations ([3] IS by numerical integration of Eq.lO forward in 
time on a rectangular grid in {E, R} space, with the density 
profile recalculated at each time step from Eq.©. We start 
from simple power-law distribution function 

f{E, R) = foR-^{E/Eof , c^iAGMbh) > \E\ > \Eo\ 

(9) 

Here the lower boundary Eq defines the energy corresponding 
to the radius of black hole's influence, where the gravitational 
potential becomes dominated by surrounding stars or dark 
matter rather than black hole itself. The higher boundary is de- 
termined by black hole horizon. The range of energies under 
consideration is large, say, 10^, to avoid boundary effects. The 
coiTesponding initial density profile is pj„ ~ K r^'^'^/2+p) 

We have set p = 3/4 (p cx r^^/"*), as it corresponds to a 
NFW halo lid adiabaticaly compressed near a black hole fit]. 
We scale density and radius to the values at the black hole 
influence radius rh, and consider r ^ r/j. The evolution was 
calculated for f3 — —2, i.e. dominance of circular orbits, and 
for (3 = 0, the case of isotropy. 

Results are shown on Fig. [T]for time ti taken so that an- 
nihilation plateau density pa equals 10^ in scaled dimension- 
less units [<i is related to pa by Eq.(|2|l]. The density profile 
evolves almost self-similarly if we scale Va and p simultane- 
ously to remain on initial profile curve, so the value ti may be 
taken quite arbitrary. 
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FIG. 2: _R-averaged distribution function: dotted line - circular or- 
bits (P = -co), solid line - P = —2, dotted-dashed line ~ /3 = 0. 
In the latter case f{E) is exponentially small for high \E\, while 
in other cases it tends to corresponding to constant-density 

core. Initial distribution function f{E) oc is shown by long- 

dashed line. The values are normalized to Eo ~ energy at the black 
hole influence radius. Energy Ea corresponding to current annihila- 
tion radius is denoted by vertical line. 



The results agree well with our preliminary suggestions. 
For the case f3 < p + 1 {(3 = —2 in our calculations) a 
constant-density plateau develops (Fig.[Tl solid) with density 
Pa and radius ra = [K / pa)^/^^''^+P\ 

In the opposite case a weak cusp p oc r""'' ,ji — (1/2 + /?) 
develops, which extends up to radius ra and smoothly joins 
the initial density profile (Fig.II] dot-dashed). The i?-averaged 
distribution function in the latter case is greatly depleted in 
the region \E\ > \Ea\ = GMbh/ra compared to the case of 
density plateau (Fig.|2|l. So the cusp is indeed formed by par- 
ticles with rather low energies (high apocentre radii) and high 
eccentricities: if the dominance of radial orbits over circular 
ones is large enough, then the fraction of particles on radial 
orbits is sufficient to determine the inner density profile. One 
can see that the i?-dependence of distribution function in the 
case < 7i < 1 is biased towards radial orbits in the region 
\E\ > \Ea\ (Fig.O. 

We note that the case of isotropic (and even radially 
anisotropic) velocity distribution is much more relevant to 
cosmological dark matter halos than the case of tangential 
anisotropy. The parameter /3 in the centers of simulated halos 
is about zero or slightly positive jBl and in most analytical 
models the situation is the same |la,|9[]. Indeed a relation be- 
tween 7 and (3 proposed in ifioll suggests /3 >= —0.15 for all 
realistic 7 >= 0, which would result in a cusp. However, if a 
binary black hole was present in the center of a halo, it would 
destroy the cusp and generate a core with tangential velocity 
anisotropy |4]. In this case, however, the annihilation plays 
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FIG. 3: i?-dependent distribution function for certain values of E: 
long-dashed and solid lines are for E = 3iJa, t = and t = ti; dot- 
ted line is for E = lOEa, ti. Light (oblique) lines are for /3 = —2, 
heavy (horizontal) lines are for P — 0. As the values are given for 
E > Ea, the 7?-profile has been changed by annihilation; however, 
for P — —2 (with density plateau), the annihilation rate is quite in- 
dependent of R for fixed E, so the evolution is in the magnitude, not 
the shape of the distribution function. For P — 0, one can see that 
circular orbits (with i? ~ 1) are depleted more rapidly than radial 
(7? ~ 0). 

almost no role because density in the core falls far below pa 
(the core radius is of order the binary separation radius, which 
is much greater than typical annihilation radius). 

IV. IMPLICATIONS FOR DARK MATTER SEARCH 

The annihilation flux <I> from the direction of the black hole 
usually is represented as a product of two quantities, the first 
of them depending on particle physics and the second, called 
"astrophysical factor" J, is related to dark matter spatial den- 
sity lUl: 



j = - — ^ 2nede dip^{^yp + {RQ9y). 

-n-o Pq Jo J -00 

(10) 

Here Qmax is flie detector angular resolution (the point-spread 
function is assumed to be of Heaviside form for simplicity), 
Rq is the distance from Sun to the black hole, is the dark 
matter density near the Sun. The outer integral represents av- 
eraging over telescope's angular resolution, the inner stands 
for the line-of-sight integration. 

We may rewrite this expression in terms of the "vicinity" of 
the black hole and the "background" from outside this vicin- 
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ity: 



J = 



Pq Jo 

Ro^max ■ 



Anr^ dr p{rY + Jbkg , 



(11) 



Now we note that for present-day observational capabilities 
we cannot hope to tell cusp from core, and even to resolve 
the black hole radius of influence. For example, GLAST will 
have angular resolution flmax of order 0.1° il2i] . which cor- 
responds to spatial distance of 15 pc at the Galactic center 
(Rq = 8.5 kpc). The radius of black hole influence in the 
center of our Galaxy is about 2 pc, and the annihilation radius 
is even smaller, of order 10^"^ pc |fl3|. 

Therefore the integral in ( fTTT i is split into three terms: r < 
Va is the annihilation plateau or inner cusp, Va < r < rh 
is the black hole domain of influence, and < r < Rmax 
is the rest. At each of these intervals the density is roughly 
power-law with index 7^. It is evident that if 7/^ < 3/2, then 
the most part of the integral comes from outer boundary of 
the corresponding region, and if 7fe > 3/2, the integral is 
determined by inner boundary. 

The density in the case of plateau is given by Eq.(|2]l, and 
in the case of weak cusp is well approximated by a similar 
expression: 



Pyr) ~ ■ , P^ = Pa[t) — 

Po + Pi \ r. 



(12) 



7j is the inner cusp power-law index, and is less than 3/2. 
On the other hand, the outer power-law index of the spike is 
always greater than 3/2 (in our case it equals 9/4). There- 
fore, the flux from the whole black hole domain of influence 
is determined by ra only. It appears that the form of tran- 
sition is significant: if we replace (fT2l i by a simpler formula 
p{r) — min(po,Pi)' '^^ overestimate flux almost twice. But 
the difference between models with /3 = — 2 (plateau) and 
/3 = (r^^/^ cusp) is less than 10%, which renders the pres- 
ence of weak cusp almost undetectable. Furthermore, the ad- 
dition of annihilation outside r^, as well as the "background", 
makes the difference even smaller Nevertheless, we want to 



stress that if the power-law index everywhere outside rh is 
greater than 3/2, then this area contributes little to the total an- 
nihilation flux. It is hkely to be the case even for NFW initial 
density profile, because it should be adiabatically compressed 
by baryons during the formation of the Galaxy, and the profile 
becomes steeper than r~^/'^ Hi' 



V. CONCLUSION 

We have reconsidered the problem of dark matter annihi- 
lation around a black hole, for the case of arbitrary velocity 
anisotropy of dark matter particles. In the case of circular 
orbits the result is well-known: a constant density core of ra- 
dius Ta and density pa develops, which smoothly joins the 
initial profile [Eq.(|2|]. However, if the fraction of radially 
biased orbits is large enough, that is, if the anisotropy coef- 
ficient (3 is greater than —1/2, a weak cusp is formed inside 
Ta'. p ~ Pa{'i'/ra)~''^^^^^^ ■ The cusp consists of particles 
which spend most part of their orbital period outside Va, but 
since their orbits are elongated and their fraction is large, they 
contribute enough to the density indide Va- The particles with 
apocentre radii within Va are annihilated much faster in the 
latter case, especially on orbits with low eccentricities. 

However, unless we have a telescope which can resolve the 
radius Va, we cannot practically distinguish between a plateau 
and a weak cusp. We note that other dynamical processes, 
such as scattering of dark matter particles off stars, signifi- 



< 



rh 



cantly affect dark matter density at radii r 
(They tend to decrease the density if it was sufficiently steep 
initially). So the effect considered in this paper does not seem 
to play an important role in the evaluation of dark matter anni- 
hilation signal from the vicinity of black holes. It is mostly of 
terminological significance: the term "annihilation plateau" in 
most cases is misleading, and should be replaced by a "weak 
cusp". 
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